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ABSTRACT 

We present the first results from the CHANG-ES survey, a new survey of 35 edge-on galaxies 
to search for both in-disk as well as extra-planar radio continuum emission. CHANG-ES is 
exploiting the new wide-band, multi-channel capabilities of the Karl G. Jansky Very Large Array 
(i.e. the Expanded Very Large Array, or EVLA) with observations in two bands centered at 
1.5 and 6 GHz in a variety of array configurations with full polarization. The motivation and 
science case for the survey are presented in a companion paper (Paper I). These first results are 
based on C-array test observations in both observing bands of the well-known radio halo galaxy, 
NGC 4631. In this paper, we outline the observations and the data reduction steps that are 
required for wide-band calibration and mapping of EVLA data, including polarization. 

With modest on-source observing times (30 minutes at 1.5 GHz and 75 minutes at 6 GHz for 
the test data) we have achieved best rms noise levels of 22 and 3.5 /iJy beam~^ at 1.5 GHz and 6 
GHz, respectively. New disk-halo features have been detected, among them two at 1.5 GHz that 
appear as loops in projection. 

We present the first 1.5 GHz spectral index map of NGC 4631 to be formed from a single 
wide-band observation in a single array configuration. This map represents tangent slopes to 
the intensities within the band centered at 1.5 GHz, rather than fits across widely separated 
frequencies as has been done in the past and is also the highest spatial resolution spectral index 
map yet presented for this galaxy. The average spectral index in the disk is ai^^cHz = —0.84 ± 
0.05 indicating that the emission is largely non-thermal, but a small global thermal contribution 
is sufficient to explain a positive curvature term in the spectral index over the band. Two specific 
star forming regions have spectral indices that are consistent with thermal emission. 

Polarization results (uncorrected for internal Faraday rotation) are consistent with previous 
observations and also reveal some new features. On broad scales, we find strong support for the 
notion that magnetic fields constrain the X-ray emitting hot gas. 

Subject headings: ISM: bubbles - (ISM:) cosmic rays - ISM: magnetic fields - galaxies: individual 
(NGC 4631) - galaxies: magnetic fields - radio continuum: galaxies 
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1. Introduction 

We present results of radio continuum observa- 
tions of the edge-on galaxy, NGC 4631, using the 
Karl G. Jansky Very Large Array (hereafter, the 
Expanded Very Large Array, or EVLA) in its C 
configuration. These are the first results from a 
new survey. Continuum Halos in Nearby Galaxies 
- an EVLA Survey (CHANG-ES), whose motiva- 
tion, science goals, galaxy sample criteria, and ex- 
pectations for improvements over pre vious surveys 



are d escribed in a companion paper (jlrwin et al 



2QI2L hereafter. Paper I). The sample consists of 
35 nearby edge-on galaxies and observations are 
being carried out in two frequency bands (1.5 GHz 
and 6 GHz, i.e. in L-band and C-band, respec- 
tively) in all polarization products. The survey 
includes observations in three EVLA array config- 
urations (B, C, and D) but the test data presented 
here were obtained at both frequencies in the C 
configuration only. 

In this second paper, our goals are to outline the 
observations and data reduction procedures used 
for the test data (Sect. [2]) and to present the initial 
scientific results for NGC 4631 (Sect.©. We pay 
particular attention to the data reduction proce- 
dures, especially with respect to differences that 
are introduced by using the wide EVLA frequency 
bands that are now available. 

The galaxy, NGC 4631, was chosen because of 
its extensive radio continuum halo which has been 



known for some time (|Ekers fc Sancisi II1977I ) . For 
a recent su mmary of previous observations of this 



galaxy, see llrwin et al. I (|2Qllf ). In this paper, we 
refer to the region, 0.2 ^ z ^ 1 kpc, as the disk- 
halo interface and use 'halo' for emission on larger 
scales {z ^ 1 kpc). 'High-latitude' or 'extrapla- 
nar' are also used to describe either of these com- 
ponents. It is worth keeping in mind that halos are 
not necessarily smooth since substructure is gener- 
ally observed, depending on the spatial scales that 
are probed. 



Queen's University, Kingston, ON, Canada, K7L 2T3, 
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2. EVLA Observations of NGC 4631 
2.1. Observations 

Observations of NGC 4631 were carried out on 
11 Nov. 2010 with the pointing center set to the 
center of the galaxy at RA = 12^ 42"^ 08^01, DEC 
= 32° 32' 29^/4 The observing set-up is summa- 
rized in Table [TJ All EVLA observations (includ- 
ing continuum) are now done in spectral line mode 
which facilitates the excision of radio frequency in- 
terference (RFI) from the data, effectively enables 
broader uv coverage via multi-frequency synthesis 
techniques (Sect. 12. 4p and permits the extraction 
of spectral index information (Sect. E^jlEl- 

At the time of the observations, the WIDAR 
correlator was restricted to a maximum of 8 spec- 
tral windows (spws) per base-band, i.e. per AC 
or BD intermediate frequency (IF). Each spec- 
tral window has its own bandpass response, and 
contains 64 channels. With AC and BD IFs con- 
tiguous in frequency, the result was 1024 spectral 
channels centered at 1.5 GHz (L-band) for a total 
bandpass of 512 MHz, and 1024 channels centered 
at 6.0 GHz (C-band) for a total bandpass of 2.048 
GHz. 

Clearly, at the lower frequency, the achievable 
bandwidth cannot be higher than the frequency 
itself, but the L-band bandwidth (more so than 
C-band) also suffers from RFjl, restricting the 
bandwidth further. In fact, although we avoided 
known RFI at the ends of the band, we found 
that the L-band RFI, especially in the 1.52 to 
1.64 GHz region within our band, to be persis- 
tent throughout our observations. Together with 
other flagging, this reduced the effective L-band 
bandwidth to approximately 300 MHz and has 
prompted a modification of the correlator set-up 
for subsequent observations. The strength of the 
RFI in L-band achieved a maximum value which 
was 300 times that of 3C286 (the flux calibrator) 
so that even the sidelobes of the RFI swamped 
that of our astronomical sources at some frequen- 
cies. Experimentation with flagged and unflagged 
data showed that it was not, in general, possible 
to salvage the long-spacing uv data while flagging 



^It will also allow an analysis of Faraday rotation, to be 
described in future papers; see Paper I for related science. 

list of known RFI is given at 
https : / / science . nrao . edu/ facilities/ evla/ observing/ rf i. 



only the short-spacings; once RFI was detected, 
all uv data had to be removed over the time and 
frequency ranges that were affected. RFI at C- 
band was also present, but much less severe than 
in L-band. 

Observations were carried out in the standard 
fashion, including scans of a flux calibrator (which 
was also used as a bandpass calibrator as well as 
for the determination of absolute polarization po- 
sition angle), a 'zero-polarization' calibrator (to 
determine the instrumental polarization) and a 
phase calibrator. The flux calibrator was observed 
in 3 independent scans at each frequency for these 
test observations. Only a single scan is neces- 
sary, in principle, but these data were intentionally 
'over-calibrated' for consistency checks. The zero- 
polarization calibrator was observed once at each 
frequency. 

The L-band observation of NGC 4631 was car- 
ried out in a single time block containing one on- 
source scan flanked by phase calibrator scans (plus 
the flux and zero-polarization calibrator). The C- 
band observation of NGC 4631 with its calibrators 
was similarly observed in a single time block, but 
separated into 3 individual scans interspersed with 
phase calibrator scans. Observations at each fre- 
quency were carried out within single blocks of 
total observing time in order to minimize time 
spent on overheads. For the regular data, how- 
ever, attempts are being made, where possible, 
to split the observing time for given galaxy at 
each frequency into two observing blocks that are 
widely separated in hour angle; more time is then 
spent on overheads but the result should be im- 
proved uv coverage with a better-behaved beam, 
as well as lower off-axis instrumental polarization 
(Sect. ESI). 

2.2. Total Intensity Calibration 

The data were reduced using the Common As- 
tronomy Software Applications (CAS A) packag^, 
unless otherwise indicated. As improvements to 
CASA are frequently being made during EVLA 
commissioning, a number of versions of CASA 
were employed in the data reductions, in particu- 
lar versions that were extant over the period May 
2011 to July 2011. The following steps were car- 
ried out for both the L-band and C-band data sep- 



Available at |http : //casa . nrao . edu] 
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arately. 

On-line Manning smoothing is no longer an op- 
tion at the EVLA, so all data were first Hanning 
smoothed to minimize Gibbs ringing in frequency, 
a step which significantly ameliorated RFI detec- 
tion throughout the band. As a result, the fre- 
quency resolution is twice the channel width, the 
latter listed in Table [fl. A significant amount of 
time was then spent fiagging RFI. Experimenta- 
tion with automatic fiagging routines had not yet 
provided results that were as dependable as man- 
ual fiagging, requiring that fiagging be done by 
hand; however, such routines are improving and 
should shorten the data reduction process in the 
future. 

The fiux density, I^y, (Stokes I) of the fiux cali- 
brator, 3C 286, was first set for the frequency, 
of each spectral window using a known model and 
the 'Per ley-Taylor- 99' calibration scal^. Since 
each spectral window of each antenna has its own 
bandpass response, bandpass calibrations were re- 
quired, again using 3C 286. However, since the 
data are vector- averaged when producing an aver- 
age bandpass for each spectral window, it is nec- 
essary to first carry out an initial gain and phase 
calibration with time; applying the solution en- 
sures that phase variations that may occur over 
the different scans of 3C 286 when it is observed 
through different elevations do not result in decor- 
relation when vector- averaging occurs. With the 
initial gain/phase solutions applied, bandpasses 
were then determined, one for each spectral win- 
dow of each antenna, and the initial gain and 
phase solution is no longer used. 

The gain and phase were then calibrated with 
time for each of the fiux, zero-polarization, and 
phase calibrators, where an antenna-based solu- 
tion is found for each spectral window and each 
separate scan, followed by bootstrapping of the 
fiux from the fiux calibrator to both the zero- 
polarization (OQ 208) and phase (J1221+2813) 
calibrators for each spectral window. At L-band, 
the fiux densities of OQ 208 and J1221+2813 var- 
ied systematically over the band from 0.72 to 1.16 



^If corrections are required for antenna position or antenna- 
based delays, they would be carried out prior to the Han- 
ning smoothing step. 



Jy (positive spectral inde:x[l with an error of about 
0.3% on each spectral window value) and from 0.55 
to 0.53 Jy (negative spectral index with an error 
of about 0.7% on each spectral window value), re- 
spectively. At C-band, the fiux densities of the 
two calibrators (in the same order) varied from 
2.4 to 1.9 Jy (negative spectral index with an er- 
ror of about 0.4% on each spectral window value) 
and from 0.46 to 0.44 Jy (negative spectral in- 
dex with an error of about 0.5% on each spec- 
tral window value), respectively. The changing 
spectral index of OQ 208 (positive in L-band and 
negative in C-band) is in good agreement with 
the known spectrum which peaks a round 4 GHz 



(^Cenacchi. Dallacasa. fc Orfei II2006L or the Nasa 



Extragalactic Database). 

The bandpass and gain/phase solutions as a 
function of time were then applied to all sources, 
antenna by antenna, and spectral window by spec- 
tral window, with the gain/phase solutions from 
the phase calibrator interpolated in time when ap- 
plied to NGC 46310. 

2.3. Polarization Calibration 

Polarized intensity, P, and the polarization an- 
gle on the sky, Xi require calibration of Stokes Q 
and U. To initially set Stokes Q and U for 3C 286, 
the absolute position angle on the sky is assumed 
to be constant at x = 33 deg, but the percentage 
polarization varies with frequency. We obtained 
detailed EVLA data on this calibrator (R. Perley, 
private communication) and fit a curve through 
the L-band data of the form, 



= -1.443 X 10 



-6 ,,2 



0.006026 + 3.858 (1) 



where Pj^ is the percent polarization at frequency, 
u (MHz). In C-band, only a linear fit was neces- 
sary, 

= 0.0001583 10.57 (2) 

{u in MHz). The estimated error, including mea- 
surement errors and fitted curves is of order 0.3%. 



See |http : //www . via . nrao . edu/astro/calib/manual/baars . html [ 



^In practice, application of the cahbration tables revealed 
that more flagging for RFI was required, necessitating a 
second iteration through the calibration. The fluxes quoted 
in this section apply to the flnal iteration. At various stages 
in the data reduction, occasionally additional flagging of 
RFI was carried out. Wherever such flagging affected the 
calibration or imaging, the necessary steps were redone. 
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Q and U were then calculated at the frequency of 
each spectral window via 

= (P./lOO) * I, * cos{2x) (3) 
= (P,/100)*I,*5m(2x) (4) 

where the total intensity, I^y, has previously been 
set via a known model (Sect. 12. 2p . 

There is leakage between the right circularly po- 
larized (R) and left circularly polarized (L) feeds 
at the EVLA and this leakage must be deter- 
mined in order to calibrate the cross-terms (RL, 
LR) from which Q and U are calculatecjl. The 
leakage terms (referred to as the 'D' terms) were 
determined from the 'zero-polarization' source, 
OQ 20dl after first applying the previously deter- 
mined gain/phase and bandpass solutions. 

Finally, the absolute polarization angle on the 
sky (the R-L phase difference) was calibrated from 
the known position angle of 3C 286, after first ap- 
plying solutions for gain/phase, bandpass, and D 
terms. 

All corrections were then applied to the cali- 
brators and the source (gain/phase, bandpass, D 
terms, and absolute position angle corrections). 

The polarization calibration described in this 
section corrects for ionospheric Faraday rotation, 
provided that this rotation is constant with time 
and position in the sky. As neither may be 
the case, it is important to determine the mag- 
nitude of such variations. At present, there is 
no CASA routine to do this, but there is (the 
task, TECOR) in the Astronomical Image Pro- 
cessing System (AIPS). We therefore copied a 
small amount of L-band data from each calibra- 
tor, ensuring that all observing times were repre- 
sented in the copied data set, into AIPS for this 
purpose. (Corrections for the source are assumed 
to be the same as those for the phase calibrator, 
given their close proximity in the sky.) TECOR 
derives corrections for ionospheric Faraday rota- 
tion from maps of total electron content obtained 
from NASA's crustal dynamics data interchange 
system (CCDidB which cover the time range of 
the observations. These measurements provide a 
crude correction since the spatial resolution is typ- 
ically coarser than the source size. Therefore, if 



significant fluctuations exist on smaller scales, we 
cannot correct for them. 

In L-band, we find that the differential iono- 
spheric Faraday rotation ranges from a minimum 
of 2.365 to a maximum of 2.805 rad m~^ over all 
sources, positions, and times. This result corre- 
sponds to a variation of 0.44 rad m~^, or 0.74 deg 
at the high frequency end of the band and 1.45 deg 
at the low frequency end. Thus a typical error 
introduced by ignoring time variable or position 
variable ionospheric Faraday rotation is about 1 
deg and we therefore have made no additional cor- 
rection for Faraday rotation at L-band beyond the 
non-differential term that has already been taken 
into account in the normal data reduction. More- 
over, since Faraday rotation cx and since the L- 
band and C-band observations were close together 
in time, no additional correction was necessary in 
C-band either. 

Finally, the accuracy of wide field linear polar- 
ization imaging is limited by variations in instru- 
mental polarization as a function of angle from 
the pointing center; that is, the corrections that 
have now been applied are accurate only for the 
on-axis position. The off-axis instrumental polar- 
ization for a 'snapshot' observatiorJ^ is estimated 
to be less than a few percent within the half- 
power point of the primary beanj^ and increases 
at larger angles. In this paper, our polarization re- 
sults (Sect. 13. 2p are shown only over regions that 
are within the primary beam full width at half 
maximum (FWHM). In addition, when sources 
are tracked, the source parallactic angle rotates 
through the primary beam, minimizing the off- 
axis polarization errors. Hence for these data, this 
effect is negligible. Nevertheless, for future wide- 
field polarization mapping, a correction will likely 
be required; a CASA algorithm is currently being 
developed to aid in this process. 

2.4. Imaging 

Our goal is to form single I, Q, and U maps 
of NGC 4631 from each of the L-band and C- 
band data; maps of linear polarization intensity. 



8Q = (RL + LR)/2, U = (RL - LR)/2i. 

^ http : / / evlaguides . nrao . edu/ index . php?t itle=Category : 
^ ^http : //cddis . gsf c . nasa . gov/] 



■'^■'^We take this to mean any observation short enough that 

Earth rotation does not fill in uv tracks. 
■"^■^See the Observational Status Summary at 



http : //evlaguides . nrao . edu or earlier references such as 



Polarim^ ^g j . //yyy . yia . nrao . edu/astro/guides/vlas/current/node35 . html | 



and lBriskenI (|20Q3l V 
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P, and magnetic field position angle, can then 
be formed from Q and U. As has always been done, 
the uv data are first Fourier Transformed (FTed) 
into the map plane which results in an image (the 
'dirty image') which is the convolution of the FT 
of the antenna distribution in the uv plane (the 
'dirty beam') with the source sky brightness distri- 
bution. The dirty beam is then deconvolved from 
the dirty image (i.e. it is 'cleaned') to f ind an 



imag e of the source. The Clark method ([Clark 



198Q[ ) of separating the cleaning process into ma- 
jor and minor cycles, treating each Stokes plane 
separately, was used. 

The broad bandwidths of the EVLA WIDAR 
correlator, however, introduce some new imaging 
requirements and considerations that were, for the 
most part, unnecessary with the VLA. These are 
described briefly here: 

a) Bandwidth smearing: If the observations had 
not been in spectral line mode, then the broad fre- 
quency band would have introduced severe chro- 
matic aberration into the image; this aberration 
worsens with distance from the map center. How- 
ever, since spectral line mode is being universally 
employed at the EVLA, bandwidth smearing is 
instead restricted to the spectral resolution of the 
data (by using the multifrequency synthesis tech- 
nique, see next section). For example, for our 
data at the center of either band and at a distance 
from the map center that is twice the FWHM of 
the primary beam, the loss in amplitude due to 
bandwidth smearing is only about 1.5%, or essen- 
tially negligible where most of the emission occurs. 

b) Multi-frequency synthesis (mfs): The posi- 
tions of the antennas in the uv plane are mea- 
sured in wavelengths and those positions will vary 
as the frequency changes across a given band. 
Consequently, if each spectral channel is grid- 
ded separately, there is an improvement in uv 
coverage over a case in which channels are av- 
eraged together. The clean algorithm employed 
in CASA utilizes this multi- frequency synthesis 
technique (jConwav. CornwelL fc WilkinsonTlllQQOl 
also called 'bandwidth synthesis'). 

c) Widefield imaging: The sensitivity of the 
array to emission that is far from the map 
center and even outside of the telescope's pri- 
mary beam means that the array itself is non- 
coplanar to such emission. Thus, the 'w term' 
in the FT (the coordinate in the line-of-sight 



direction towards the source), which can be ig- 
nored near the map center, becomes important. 
CASA clean employs the 'w-projection algorithm' 
(|CornwelL Golap. fc Bhatnagar 1 12008 ) to correct 
for this effect. 

d) Spectral index fitting: With wide bands, the 
spectral index of the source itself, which is spa- 
tially variable, need s to be fitted during imaging 
(|Sault fc Wieringa I [1994). In CASA's clean algo- 
rithm, this is accomplished by assuming that the 
emission at any position can be described by. 



I. 



I.. ( ^ 

.^0 



(5) 



where is a reference frequency within the band 
(here taken to be the band center), ly^ is the spe- 
cific intensity at the reference frequency, a is the 
spectral index and j3 describes the spectral cur- 
vature. The function is expanded in a Taylor se- 
ries about resulting in a polynomial fit to the 
spectrum. The fit includes all unflagged channels 
across any given band and is not done on a 'per 
spectral window' basis. In addition, the fitting is 
done using all data so that the resulting parame- 
ters correspond to a single spatial resolution (for 
example, the a map is not formed such that the 
spatial resolution would vary across the band it . 
lR,au fc CornW n ("2011") outline the implementa- 
tion of spectral fitting and provide examples as to 
how the image can improve with increasing num- 
bers of Taylor terms, ut- 

e) Multi-scale clean: Improvements have been 
made over the traditional clean which modeled 
all emission as a collection of point sources. The 
clean implementa tion in CASA n ow allows for a 
multi-scale clean (|Cornwell Il2008[ ) which assumes 
that the emission can be modeled as a collection of 
components over a variety of spatial scales. These 
scales are typically chosen to span the range from 
an angular size of 'zero' (corresponding to a point 
source, or the classic clean analogy) to a scale of 
order the maximum spatial scale expected to be 
present in the emission. Since the latter is not al- 
ways known in advance, in practice, the maximum 
scale size is increased systematically until no more 
flux is actually cleaned out in that largest scalJ^. 



^^At the time of writing, this is the only CASA-based fitting 

function that is available. 
-•^^If the fiux becomes negative in the largest scale, the trial is 
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For details of the CASA implementation of the 
combined multi-sc ale, multi- frequen c y syn thesis 
(ms-mfs) clean, see i Rau fc Cornwell 
f) Primary beam: The primary beam (PB) varies 
from one end of the band to the other and, as a 
result, nowhere in the field of view is there a 'null' 
as the re would be for a mono chromatic primary 
beam (|Bhatnagar e t al. ]l2Qll[ ). Consequently, it 
is the frequency- averaged PB (PB^^^) that the fi- 
nal cleaned image must be corrected for. More- 
over, since the frequency dependence of the PB 
imposes its own spectral index across the field of 
view, all spectral index maps (Sect. 12. 6p must be 
corrected for this effect as well. 

NGC 4631 was mapped taking into account all 
of the above points. 

At both bands, a very wide field was first 
mapped (3 degrees for L-band, or 6 times the pri- 
mary beam FWHM; 34.1 arcmin for C-band or 

4.6 times the primary beam) to determine how 
large an image needed to be mapped to include 
all sources of significant flux in the clean. Smaller 
fields were then mapped completely, as required. 
For the multi-scale clean, we used typically 4 to 6 
equally-spaced spatial scales, the first one always 
corresponding to point sources as in the classic 
clean. For example, at L-band, scales up to 4.2 
arcmin were employed, where the largest scale was 
chosen according to the considerations of point 
e) above. We also self-calibrated each data set 
once, finding that a single amplitude + phase self- 
calibration resulted in a small improvement in the 
rms map nois^H- 

After extensive experimentation, we adopted 3 
sets of parameters representing 3 different spatial 
resolutions for each frequency and display them, 
along with some ancillary data, in Figs. [1] and 
[2j the corresponding map parameters are given in 
Table El 

The theoretical noise values for these observa- 
tions (including confusion estimates) are 18 and 

3.7 /iJy beam~^ at L-band and C-band, respec- 
tively. Our lowest rms noise values are only 1.6 
and 1.3 times higher for the total intensity im- 
ages. There are, however, several residual linear 
cleaning artifacts in the L-band images that are 



visible immediately to the east and west of the 
galaxy's disk (see Fig. [Th in particular). The con- 
clusions that we present in this paper are not af- 
fected by these artifacts. The disconnected emis- 
sion seen approximately 2 arcmin above and below 
the plane in Fig. [2t is likely real and part of the 
halo; missing low-order spacings contribute to the 
well-known appearance of the galaxy's disk emis- 
sion sitting on a negative bowl. 

Note that, as is usually done, all displayed to- 
tal intensity images are not corrected for the pri- 
mary beam in order to show images with uni- 
form noise. When numerical results are given, 
however, the results have been corrected for the 
frequency- averaged primary beam (PB^^^). Note 
also, that we carried out some subsequent analysis 
using AIPS. 

2.5. Polarization Imaging 

Stokes Q and U maps were formed with the 
same sets of parameters as the total intensity im- 
ages, except that fewer spatial scales were required 
during cleaning, given the weaker emission in these 
maps. The relevant parameters for these images 
are given in Table O In principle, the rms noise 
on the Q and U maps should be the same as the 
I map; however, the lower signal in the cross- 
hands means that residual errors associated with 
cleaning the high intensity emission are no longer 
present. At L-band, our lowest rms noise is 1.2 
times the theoretical value and at C-band, our 
lowest noise matches the theoretical value. Note 
that, in order to obtain sufficient signal-to- noise 
(S/N) in the polarized emission, all spectral win- 
dows were used. 

We formed polarized intensity images, P= 
a/Q^ + U2, for each of the uv weightings of Figs.[T] 
and [21 correcting for the bias introduced by the 
fact that P images do not ob ey Gaussian statistics 
(Simmons & Stewart 1 Il985|; [Vaillancourt 20061). 
Maps of the polarization angle of the electric field 
vector, X = (l/2)arctan(U/Q), were also formed 
forpoints for which the linear polarization was > 
4cr% When displayed, we have rotated x by 90 
degrees to illustrate the observed (uncorrected for 
internal Faraday rotation) direction of the mag- 
netic field. We show the results for only a single 



aborted and the maximum scale is decreased until positive 
flux is cleaned in all scales. 
-•^^A phase-only self-calibration did not improve the results. 



-"^^The cutoff for x was relaxed slightly from 5cr in order to 
show the trends in vector orientation with position. 
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representative data set in the two frequency bands 
(Robust = 2 for L-band and Robust = + uv ta- 
per for C-band) in Fig. [jPl . Consistent with the 
total intensity images (Sect. 12. 4p the displayed 
maps are shown without primary beam correc- 
tions. 

We then carried out a PB correction using the 
same frequency- averaged primary beam, FBavg, 
as was used for the total intensity images. This 
correction is good to zeroth order, but does not 
take into account first order effects which would 
result from differences between the right and left 
circularly polarized PBs. First order corrections 
are not yet implemented in CASA but we expect 
them to be small for this data set. For example, if 
the difference between the right and left circularly 
polarized PB voltage response is of order 10% at 
a distance from the map center corresponding to 
30% of the primary beam (17.4 arcmin and 4.6 
arcmin at L-band and C-band, respectively), then 
the error in Q or U would be 3% of the Q or U 
intensity at that point. 

PB corrections of the linearly polarized inten- 
sity maps shown in Fig. [3] of all emission greater 
than 5<j resulted in increases in flux of 5% and 
40% at L-band and C-band, respectively. Maps of 
X (shown in Fig. [3]) do not require PB correction. 

We then formed maps of percentage polariza- 
tion for the two bands (these also do not require 
PB correction). The maps were formed for emis- 
sion which was greater than 5cr in both linear po- 
larization intensity as well as total intensity, ensur- 
ing that the spatial resolutions are matched. We 
do not display the L-band maps since the result 
had sufficient S/N in only 20 beams and these po- 
sitions were non-contiguous spatially. The C-band 
percentage polarization map is shown in Fig. IH 
We calculated the average percentage polarization 
over the region shown by measuring the total po- 
larized flux (rather than per pixel values) and di- 
viding by the total flux, so that high values around 
the periphery of the emission do not skew the re- 
sult; we find an average of 7.0% for the result cor- 
responding to the region displayed. This region is 
mainly in the disk where some Faraday depolar- 
ization may be occurring (Sect. 13. 2p . 



2.6. Formation of Spectral Index Maps 

As indicated in Sect. 12.41 point d)^ the wide 
bands that are now available at the EVLA make it 
possible, in principle, to create spatially resolved 
spectral index maps, a, (and curvature maps, /3, 
if there were sufficient signal-to-noise) over a sin- 
gle band with a single observation and common 
spatial resolution. One could form a map of a 
for L-band and for C-band individually, and then 
form such a map between the two bands, assuming 
a constant spectral index between the two bands. 

For observations at two frequencies with a sin- 
gle EVLA array, as for our C- array test observa- 
tions, however, there are some limitations. First, 
as specified in Paper I, we do not detect spatial 
scales at L-band and C-band that are greater than 
16.2 and 4 arcmin, respectivel\l^. The halo of 
NGC 4631 at both bands has previously been ob- 
served to extend to a pproximately 12 arcmin (see 



for ex ample, Fig. 9 oflHummel. Beck, fc Dettmar . 
(1991) and Fig. 3 of lCoUa fc Hummel I (11994*)). It 
is therefore clear that not all flux has been de- 
tected at C-band and, as Fig. [T] illustrates, the 
broadest scale halo emission has not been detected 
at L-band either. Future D-array observations, 
with additional shorter spacings, should remedy 
this situation at L-band. 

Using previously published total flux esti- 
mates, we can normally quantify the amount of 
flux that has been missed in our observations. 
However, at C-band, a wide range of 4.85 GHz 
fluxes have been reported, likely reflecting sim- 
ilar problems with spatial scales, for example, 
S4.8GHZ = 438 mJy, 243 mJy (see NED) and 
480 mJy (|Golla fc Hummel I Il994l ). Our primary- 
beam-corrected flux in C-band is Sqghz = 180 
mJy. Using the aver a ge glo bal spectral index of 
iHummel fc Dettmar I ( 1990[ ). aavg = —0.73, our 
frequency- adjusted flux would be S^^sghz = 210 
mJy, indicating that there is significant flux miss- 
ing from the C-band maps. We thus do not show 
our C-band spectral index maps in this paper. 



■"^^Note that the CASA package did not allow for variable 
vector lengths, hence they are drawn as constants in these 
figures. 



-•^^ These values are for full-synthesis observations. For 
short observations of ^ one minute duration (snapshots) 
these values could reduce to as little as 8 and 2 ar- 
cmin, respectively (see Observational Status Summary at 
[http://evlaguides.nrao.edu). Our observations (see Ta- 
ble ^ are intermediate between a snapshot and full syn- 
thesis. Consequently, the spatial scales quoted in the text 
are likely upper limits. 
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At L-band, we do form spectral index maps. 
The primary-beam-corrected flux from these C- 
array observations is Si.^ghz = 935 mJy, or 77% 
of the VLA D-ar ray VLA fluxes of 1 22 Jy at 
u = 1.49 GHz (|Hummel Dettmar I \l9m or 
1.29 Jy from the Westerbor k Synthesis Radio Te le- 



scope dit u = 1.365 GHz (|Braun et al. ll2QQ7l i^ 



Since the missing flux is on the largest scales, we 
form spectral index maps at L-band beginning 
with the highest resolution image (Fig. [T])a, and 
then consider how this and other errors affect the 
result. Since the L-band 500 MHz bandwidth is 
much narrower than the frequency spacing which 
has been used in the past to calculate a (e.g. 610 
MHz to 1412 MHz, or 4.75 GHz to 10.7 GHz, see 
Werner 19881 ). the values of a presented here can 



be thought of as tangents to the frequency spec- 
trum at each spatial position in L-band. 

Spectral index and curvature maps were formed 
from ratios of the Taylor expansion coefficient 
maps which describe Eqn. I^Pl . after first correct- 
ing these coefficient maps for PB^^^ and then in- 
troducing a 5(7 total intensity cut-off. It has been 
shown that increasing the number of terms to the 
polynomial fit will generally in iprove the result 
subje ct to S/N considerations teau fc Cornwell" 



201 If ). For example, fitting both a and /3 allows 



a to be well-determined at the reference frequency 
even though /3 may not have sufficient S/N in each 
synthesized beam to be presented as a map. In 
Fig. [5] we show spectral index maps corresponding 
to the high resolution L-band data (Fig. [TJi) for 
a straight spectral index assumption (a) and then 
allowing for curvature within the band (b). 

The "mottled" appearance of the two maps is 
due in part to the fact that point-to-point varia- 
tions are significant. For example, after subtract- 
ing the two spectral index maps, the rms of the 
result is cr ^ 0.4 (omitting approximately a beam 
width around the perimeter of the emission) sug- 
gesting that the per-pixel error from the fitting is 
of this order. There are, however, variations in 
the map that exceed this uncertainty. In addi- 
tion, there are some global trends in a which are 
very similar between the two maps and the aver- 



i^The NRAO VLA Sky Survey (NVSS) flu x of 981.6 mJy 
has b een underestimated for this galaxy (seel Condon et ahl 
12002 ). 

= TTl/TTO and /3 = TT2/TT0 - a(a - l)/2, where 
TTO, TTl, and TT2 are maps of the first 3 Taylor terms. 



age spectral indices (a, averaged over the maps) 
are in very good agreement, differing by only 0.01. 
Clearly, the error bar on the spectral index de- 
creases with the number of points that are aver- 
aged, i.e. with increasing S/N, as one would ex- 
pect; we will return to this point in Sect. 13.31 

Regarding the effect of missing broad scale flux, 
our measured a will not be affected provided the 
broad scale emission has the same spectral index 
as the disk and contributes a constant intensity 
at each position. It is known, however, that the 
halo spectral index, which we associate with the 
broad scale emission, is steeper than that of the 
disk. A variety of measurements have been made 
over different frequency ranges, but typically, off- 
major- axis values fall in the range, a = —0.9 
—2 (see discussion in lHummel fc Dettmar |[l990[ ). 

Taking the extreme case of a = —2, we can 
determine how the addition of a 'plateau' of emis- 
sion of this spectral index to our data would af- 
fect our measured values of a. If the halo emis- 
sion occupies an emission ellipse of major x minor 
axis extent, approxim ately 15 arcmin x 12 arcmin 



(|Hummel et al. Ill988i) . then the emission region 



of Fig. [5] occupies approximately 20% of this area. 
Scahng the missing flux of Fig. [TJi to this area, as- 
suming that the broad scale flux uniformly fills the 
emission ellipse, then the missing emission corre- 
sponds to a plateau of 0.12 mJy beam~^, on aver- 
age, over the disk at the central frequency. 

Using the measured average intensity from 
Fig. [TJi over the same region as the spectral in- 
dex (taking /3 = 0) the addition of such a steep 
spectral index plateau would adjust our measure- 
ment of a from -0.84 to -0.88 (Aa = 0.04). This 
is a worst-case scenario unless there are missing 
spatial scales that are smaller thdiH 11 arcmin (the 
major axis size). 

Spectral index, a, and curvature, /3, maps for L- 
band were also formed for the remaining mid and 
low resolution maps corresponding to Fig. [Ud and 
c, respectively and are shown in Fig. [6l For these 
resolutions, we find a = —0.82 and a = —0.87, 
respectively. The curvature, averaged over the 
emission regions for these two resolutions is 1.9 in 
both cases, agreeing to within ±0.1 with the result 
at high resolution. If we consider fitting errors, 
possible errors from missing flux, and variations 
between the different uv weightings, our average 
spectral index is a = —0.84 ± 0.05. 
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As mentioned earlier, it is worth stressing that 
the wide bands used here allow for considerable 
flexibility in 'tuning' the spatial resolution of the 
array by appropriate use of uv weighting. This 
means that the spatial resolutions of L-band and 
C-band data can actually be matched exactly, 
even with observations in only a single array con- 
figuration. Therefore spectral index maps can also 
be made between L and C-bands from observa- 
tions in a single array configuration. Although 
we formed such maps for NGC 4631, we defer the 
results until all missing C-band fiux is recovered 
through additional observations (Paper I). 

3. Results and Discussion 

These new EVLA observations, even though 
they were obtained during a limited test run in 
a single array configuration (Table [1]), have re- 
vealed a number of new features not before seen in 
NGC 4631. In this section, our intention is not to 
make an exhaustive list of all features, but rather 
summarize the more obvious or new results and to 
place them in the context of previous knowledge 
of this galaxy. 

3.1. Total Intensity Images and Compari- 
son with other Wavebands 

The total intensity images displayed in Figs. [T] 
and [2] show numerous extensions away from the 
plane at both frequencies. 

At L-band, in particular, two of the more 
prominent features appear in projection as loops 
or partial loops (although they could represent 
shells) that have been outlined in red in Fig. [TJd. 
Although we have noted a few cleaning artifacts 
in this image to the east and west of the major 
axis of the galaxy (Sect. 12. 4p . the loops do not 
fall within the affected regions and they are ei- 
ther of sufficient signal-to- noise and/or persistent 
over various uv weightings (some not displayed 
here but formed during experimentation with uv 
weightings as described in Sect. 12. 4p that we con- 
sider them to be real features. 

The diameter of the larger northern loop, which 
is only partial in Fig.[T|3, is 6.3 kpc. It is also visi- 
ble in the lower resolution Fig. [ih where it appears 
to 'frame' the companion dwarf elliptical galaxy, 
NGC4627 which is 2.6 arcmin to the NW of the 
galaxy's center. This is likely a coincidence, since 



we will show below that the loop is related to the 
hot X-ray gas. 

The diameter of the smaller more easterly loop 
is 2.6 kpc and this loop becomes 'filled in' in the 
lower resolution Fig. [ih. The eastern side of this 
loop and the weaker emission that extends up- 
wards fro m it corresponds to the HI ' worm' ob- 
served bv lRand fc van der HulsTI (|l993h . The lat- 



ter feature sits above a significant HI supershell 
found by those authors whose diameter (3 kpc) is 
approximately the same as our eastern loop. It is 
likely that these features are related and are asso- 
ciated with star forming regions that are present 
in the underlying disk. The closest prominent star 
formation complex is located at RA ^12^ 42"^ 22®, 
DEC ^ 32° 32' 42''; it is marked with a cross in 
the Ra map of Fig. [St. 



NGC 4631 (Hummel et al. Ill988: Hummel & Dettmar 


1990: "Hummel. Beck 


& Dettmar 1991: Golla & Hummel 


1994: Dumke et al. Ill 


L995: Goha 1999: Braun et al. 1 


2007:lKrausell2009l 


Heald et al. l2009h does not 



reveal these loop s, alt hough the continuum image 
(|2QQ7) shows hints of the larger 



in iBraun et al. 
central loop. Since we see extraplanar structure 
on all scales in this galaxy, detecting specific fea- 
tures requires appropriately weighted uv data over 
scales that match the feature of interest. The spa- 
tial scales highlighted by our new data have now 
revealed these loops and the other arcs and fila- 
ments that we see extending from the plane. It 
is possible that some of the other broad exten- 
sions (for example, at RA ^ 12^ 42"^ 12% DEC ^ 
32° 31' 30") would resolve into loops if viewed at 
higher resolution and sensitivity. 

The close connection between X-ray emission 
an d the radio ha l o, first point ed out for this galax y 



an d tne radio na l o, nrst pomt ed out tor tnis galaxy 
by IWang et al. I (|l995h and IWang et al. I (|200lh . 

finds stronger support in these observations. In 
Fig. [71 we show the soft X-ray emission from 
Fig. [TJi, enhanced to show extra-planar emission, 
together with the low resolution L-band image 
from Fig. [Th. In this figure, we see that the X- 
ray emission forms a loop interior to the northern 
radio continuum loop, suggesting that hot gas may 
be confined by magnetic pressure. It is interesting 
that there is even a smaller vertical X-ray protru- 
sion (RA ^ 12^ 42"^ 07% DEC ^ 32° 36') above the 
gap on the north-eastern side of the radio contin- 
uum loop. A close relationship has been found be- 
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tween radio continuum and X-ray emission in th e 



outflow of NGC 253 as well (iHeesen et al. ll2Qll[ ). 



We can determine whether such confinement 
is feasible, assuming a simple case of unity fill- 
ing factors, by estimating the minimum energy 
magne tic field strength following iBeck fc K rause i 
(|2QQ5h . Assuming a proton to electron number 
density ratio of 100, a line of sight distance equal 
to the width of the radio loop as shown in Fig. [71 
an isotropic field direction, and a spectral index 
equal to the global average disk value, we find 
Bmin ^ 8 mG- The magnetic pressure is then 



Pmag — 3 X 10 



-12 3 -gy comparison, the 



thermal pressure from the hot. X-ray-emitting gas 
is Pth 6.9 X 10-^^ (Wang et al. 2001), imply- 
ing that magnetic pressure is indeed sufficient to 
confine the thermal gas. 

At C-band, where we have higher spatial reso- 
lution (Fig. [2]), we continue to see many vertical 
extensions. Again, we do not discuss every feature, 
but point out the brightest extension which is on 
the north of the disk at RA = 12^ 42"^ 09% DEC 
= 32° 33' 10'', and best seen in Fig. EJd. This ex- 
tension is directly below the west side of the large 
northern loop seen in Fig. [TJd. This feature pro- 
trudes 1 kpc above the main emission below it, 
or about 1.6 kpc from mid-plane. The detail and 
numerous extensions observed on this map reveal 
the complexity of the disk-halo interface in this 
galaxy. 

The disconnected 'patches' of emission visible 
in Fig. [2t above and below the plane appear to 
represent brighter parts of the broader scale halo 
emission which has not been detected in its en- 
tirety in our data due to missing short spacings at 
C-band (Sect. [23]). The feature at RA ^ 12^ 42"^ 
23% DEC ^ 32° 34' 30" for example, corresponds 
to the HI worm mentioned above. 

The well-known close relationship between ra- 
dio continuum. Ha, and CO emission is also evi- 
dent from Fig. [21 to the extent that 'kinks', widen- 
ings and narrowings along the major axis are cor- 
related between these components. Radio contin- 
uum and Ha peaks are nicely coincide nt in the 
wider 'bulge' region (jlrwin et al. ]|201l') about 3 
arcmin to the east of the galaxy's center within 
which the bright star formation complex men- 
tioned above is located. 



3.2. Polarization Maps 

The linear polarization maps with observed 
magnetic field vectors are shown in Fig. [3] and the 
percentage polarization with observed magnetic 
field vectors for C-band is shown in Fig. [H Recall 
that the Q and U images have lower rms than the 
total intensity images (Sect. 12. 5p and hence some 
polarized emission can be seen outside of regions 
of total intensity. 

All maps have been corrected for foreground 
ionospheric Faraday rotation, where necessary 
(Sect. [231), but no Faraday rotation correction has 
been done for the galaxy itsel£il; therefore, espe- 
cially in the disk at L-band, we expect significant 
Faraday depolarization, as indicated in Sect. 12.51 
This is evident in Fig. [3^1 which shows little po- 
larized intensity in the disk a t L-band. Previous 
L-ba nd images (see especiallv Heald et al. iliooit 
Hummel, Beck, & Dettm^J' 1991) similarly show 
little or no polarization right in the disk. The halo, 
however, should not be as strongly affected. The 
effects of Faraday rotation should also be lower 
at C-band (Fig. [Sb) but will not be negligible in 
the disk. For example, rotation measures of up 
to ± 300 rad/m^ have bee n found fo r NGC 4631 
in a larger 85 arcsec beam (|Krause Il2004f ). corre- 
sponding to 43 degrees of rotation at 6 GHz; the 
value could be higher still in our data set, given 
our smaller beam. Again, away from the plane, 
Faraday rotation should be minor, as has been 
demonstrated previously for the halo region of the 
edge-on galaxy, NGC 5775 (|Soida et al. 11201 ll cf. 
Figs. 3 and 5 of that paper). Nevertheless, we 
prefer to be cautious about our comments regard- 
ing the magnetic field orientation in the following 
discussion. 

At L-band, the most prominent halo feature is 
a highly linearly polarized spur located at RA ^ 
12^ 42"^ 20% DEC ^ 32° 34' 30". T his spur is vis- 



ible in the linearly polarized map of iHeald et al 
(I2OO9I ) and also in the lower resolution 5 GHz 



map of iGoUa fc Hummel I (|l994[ ). The spur is di- 
rectly above the eastern radio continuum loop and 
may be related to the disturbance that has caused 
that loop as well as the HI supershell below it 
(Sect. 13. ip . However, the galaxy also displays an 
V-shaped field structure on the northern side (or 
'X-shaped' globally, see above references) and this 



^We will carry out such an analysis in a future paper. 
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spur may be associated with the eastern side of 
that larger feature. 

At 6 GHz (Fig. [3b, Fig. g]), it is interesting 
that the observed magnetic field orientation within 
the disk is parallel to the disk. Previously pub- 
lished maps have shown the magnetic fields near 
the center of NGC 4631 to be perpendicular to 
the disk (for example, see the 8 GHz image of 
Krause 20041 at 84 arcsec resolution), although 
at lar ger radii the fi elds are parallel to the disk 
again (Krau se 1120091) . Our 12 arcsec observations 



and omission of the largest spatial scales are tar- 
getting different scales and may now be revealing a 
disk-parallel component near the galaxy's center, 
although Faraday-corrected maps are still required 
to confirm this. Away from the disk, the mag- 
netic field tends to become perpendicular to the 
plane, in agreement with previous observations. 
Extra-planar polarized emission both to the north- 
east and north-west of the gal axy's center agrees 



in gen eral with the image of iGolla fc Hummel 
(|l994h at 5 GHz. 



A final point of interest is the sharpness with 
which the apparent magnetic field orientation 
changes (90 degrees within a single beam) at RA 
^ 12^ 42^ 12% DEC ^ 32° 32' 30'' (see Fig. g]). 
This change occurs at the location of a bright HH 
region complex vi sible in Fig. ^ which has been 
labelled CM67 bv lGollal ^99). If such a change 
is still visible after having corrected for Faraday 
rotation, then it lends support for the notion that 
the magnetic field is swept outwards with the out- 
fiow from specific star forming regions. Note that 
there are radio continuum spurs above and below 
this location (Fig. [2t). If the change is not in- 
trinsic but is due to Faraday rotation alone, then 
either the ordered magnetic field in the line of 
sight has changed sharply at this position or the 
electron density has changed abruptly, or both. 
In either case, CM67 is an important region for 
follow-up studies with applied rotation measure 
corrections. 

3.3. Spectral Index Images 

Our L-band spectral index maps (Figs. [5] and [6j) 
are the first that have been formed from within an 
observing band for an edge-on galaxy and repre- 
sent a tangent to the spectrum at 1.5 GHz, rather 
than an average between broadly separated fre- 
quencies as has been done in the past. 



The average spectral index is o; = - 0.84 ± 
0.05 (Sect. 12. 6p and therefore the L-band maps 
are clearly dominated by non-thermal emission. 
This res ult is typical of the values found in other 

7 ) and 



galaxies (Niklas, Klein, 
the Milky Way (Reich, 



Wielebinski"T9! 
Reich & Testori 
[strong. Orlando, fc Jaffe I l201lf ). Previous 



val- 



ues for the disk of NGC 4631 coyer a rather wide 
range (see iHummel fc Dettmar |[l990[ ). for exam- 
ple, -0.68 (from 610 MHz and 10.7 GHz mea- 
surements), -0.45 (from 327 MHz and 1.49 GHz), 
-0.69 (from 327 MHz, 610 MHz, 1412 MHz, and 
10.7 GHz data), and -0.9 (from 610 MHz and 1412 
MHz). 

Also, although the global spectral index of 
NGC 4631 has been studied extensivelv (iKlein et al. 



198^ "Sukumar fc Velusamv "'1985; 'Werner "1988'; 
Hummel fc Dettmar 1990; Po hl fc Schlic keiserJ 
1990h the only previous spectral index maps that 



have been published can be found in IWerner 
y.988) whose highest resolution was 58 arcsec. 
Our maps are the highest resolution spectral in- 
dex maps yet presented for this galaxy. 

The average, global spectral curvature is posi- 
tive {P = 1.9, Sect. 12. 6() . The most likely explana- 
tion is the contribution of thermal emission from 
the disk. To see whether this is feasible, we use 



(|l999h 



the thermal fiux density estimate of lGolla 
of Ft (4.86 GHz) = 56 mJy who makes use of 
the Ha luminosity of NGC 4631 and adjusts up- 
wards to account (roughly) for dust obscuration. 
Extrapolated to 1.5 GHz with a frequency 
dependence gives Ft (1.5 GHz) = 63 mJy with a 
variation of 2 mJy across L-band. The curvature 
that we have found from Eqn. [5] corresponds to an 
increase of 2% at the upper frequency band edge in 
comparison to a non-curving spectral index. With 
a total fiux density of 935 mJy at the band center 
(Sect. 12. 6p , a curved spectral index results in an 
increase of 16 mJy at the high frequency end of 
the observing band, over the value that would be 
observed if the spectral index were constant. This 
increase is less than the estimated thermal fiux 
density at the upper end of the band, so it is clear 
that a thermal contribution is capable of account- 
ing for the observed curvature. However, thermal 
emission contributes, on average, less than 7% of 
the fiux a t L-band, ass uming that the thermal es- 
timate of lGolla I (|l999h can be applied. 

An interesting result is the detail with which we 
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have mapped the spectral index distribution in a 
single observation in a single frequency band. As 
indicated in Sect. 12. 6| the uncertainty at any point 
is approximately ±0.4, though it becomes lower 
depending on the size of the region over which 
averages are taken. For the average spectral in- 
dex over all maps, the uncertainty is smaller, i.e. 
±0.05. Consequently, in the following, we point 
out only values of a that are significant in compari- 
son to the uncertainties which have been appropri- 
ately averaged over the region of interest and for 
which there are similar trends for all maps shown 
in Figs. [5] and [HI 

For example, contrast in the spectral index is 
observed in the region of the eastern loop (see 
dashed arc in Fig. [TJd). The contrast is most 
clearly seen in Figs. [6^1 and b. The average spec- 
tral index along the sides of the loop itself (all 
maps) is a = —1.3 and this value does not differ 
significantly from the underlying disk at the loop 
'footprints'. In the regions interior to the loop and 
just exterior to it, on the other hand, the values 
are much steeper (darker regions in Fig. [6^1 and b). 
The average spectral index interior to the loop is 
a = —2.1. Both values indicate that a thermal 
contribution will be negligible at these locations 
and, as we did in Sect. 13. l( we estimate a minimum 
energy magnetic field strength using the same as- 
sumptions as in that section but with these mea- 
sured spectral indices. The result is ^rnin ^ 17 
/iG along the loopEl. 

This is not the first time that spectral index 
contrasts or steep spectral indices have been ob- 
served in galaxies that ex perience outflows. For 
example, iLee et al. ] (2001) show that the spectral 
index of radio continuum features that are located 
along HI shells in NGC 5775 is flatter than in 
the interior and that this flattening cannot be 
due to a thermal contribution. As in NGC 4631, 
the spectral index along the loop itself is not sig- 
niflcantly different from that of the underlying 
disk at its location. Moreover, although such 
steep spectral indi ces are not generally seen in 
our own Galaxy (iReich. Re ich fc Testori ' 2004*; 
Strong. Orlando^J^Taffellioil : Ghosh et al . 20lJ), 



Heesen et al. ( 201lh also flnd regions of very 



^^By comparison, a classical calculation (|Pacliolczvk Ill970l ) 
using a heavy particle to electron ratio of 40, and frequency 
limits from 10^ to 10^^ Hz, results in Bmin ~ 10 /j,G. 



steep spectral index in the galaxy, NGC 253, i.e. 
a =—2.0 ±0.2 is seen within the outflow cone 
originating from the nuclear region of that galaxy. 
NGC 4631 is known to be experiencing widely dis- 
tributed disk-related activity rather than a loca l- 
ized nuclear starburst (e.g. see llrwin et al. lHoill ). 
but the physics may be similar. For NGC 253, 
strong electron energy losses are implied, for ex- 
ample. 

We defer the point-by-point decomposition of 
the map into thermal and non-thermal compo- 
nents to a future paper; nevertheless, it is possible 
to see evidence for the dominance of thermal emis- 
sion in at least two speciflc star forming regions, 
best seen at the highest resolution spectral index 
maps (Fig. [5j). One has been pointed out earlier 
(Sect. 13. ip and is marked with a cross on the Ha 
map of Fig. [2t as well as the high resolution spec- 
tral index map of Fig. [SJd. The spectral index at 
this position is a = -0.16 ± 0.09 (averaged from 
Fig. [SJd over the star forming region) which is en- 
tirely consistent with a thermal spectrum, within 
uncertainties. Another star forming region located 
1 arcmin to the SSW of this location and visible 
as a discrete peak in Fig. [2t also shows a similarly 
flat index, consistent with thermal emission. 

4. Conclusions 

In this paper, we have presented the flrst results 
from a new survey, called CHANG-ES, to observe 
radio continuum halos in 35 edge-on, normal spi- 
ral galaxies. The galaxies are being observed in 
all polarization products in two different bands, 
1.5 GHz (L-band) and 6 GHz (C-band), and over 
3 different EVLA array conflgurations. This is the 
flrst comprehensive radio continuum survey of ha- 
los to include all polarization products. The moti- 
vation and science goals for the survey have been 
presented in Paper I. 

Our initial CHANG-ES test observations of 
NGC 4631 have been carried out at C array alone 
and therefore are sensitive to the disk, the disk- 
halo interface, and inner halo emission; we do not 
detect the faint large scale halo emission because 
of the lack of large scale sensitivity in this array 
conflguration. Our new results demonstrate that, 
even with modest integration times, new details of 
the emission in this galaxy have emerged. 

One advantage of the wide-band, multi-channel 
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EVLA capabilities is the ability to form spectral 
index maps at a common spatial resolution within 
a single observing band in a single array config- 
uration. We have formed such maps at L-band. 
Our L-band spectral index map has the highest 
resolution yet obtained for this galaxy and is the 
first in-band spectral index map published for an 
edge-on galaxy. In addition, with appropriate uv 
weighting, the EVLA wide bands have allowed us 
to match spatial resolution between L-band and 
C-band in a single array configuration. 
Our results for NGC 4631 include: 

• At both 1.5 and 6 GHz, numerous extensions 
can be seen emerging from the disk into the halo, 
many of which have not previously been observed. 

• At 1.5 GHz, two extra-planar features appear 
to form loops in projection. The larger loop is 6.3 
kpc in diameter and located on the north side of 
the disk slightly to the west of center. The smaller 
loop, 2.6 kpc in diameter, is also on the north side 
of the disk but towards the east (see Fig. [TJd). 

• The larger loop is exterior t o the s oft X-ray 

(Fig. [71). 

) juG and 
3 X 



feature observed by I Wang et al. I (|2QQll 



The minimum energy magnetic field is 
the magnetic pressure in this loop {Pmag 
10~^^ erg cm~^) is sufficient to constrain the hot 
X-ray gas {Pth = 6.9 x 10~^^ erg cm~^). 

• The smaller loop is approximately at the lo- 
cat ion of the HI worm an d supe rshell identified 
by I Rand fc van der Huls t ' (1993). The spectral 
index within and adjacent to the loop is much 
steeper {a ^ —2.1) than along the loop itself 
{a ^ -1.3). 

• At 1.5 GHz, a spur of polarized emission can 
be seen away from the plane, consistent with pre- 
vious observations; this spur is likely associated 
with the eastern loop. 

• At 6 GHz, the polarization is higher and the 
average percentage polarization is 7% over regions 
within which the S/N was high enough to be mea- 
sured, i.e. predominantly in the disk. The appar- 
ent magnetic field orientation in the plane of the 
disk (uncorrected for Faraday rotation) appears to 
be parallel to the disk. 

• There is one position in the disk, correspond- 
ing to the HII region complex, CM67, at which the 
apparent field orientation changes sharply (Fig.|4]); 
this is a target region for future in-depth analysis 
including Faraday rotation. 



• At 1.5 GHz, the average spectral index is 
a = —0.84 ± 0.05, indicating that the emission is 
predominantly no n- thermal, on average, through- 
out the disk. However, a small thermal contribu- 
tion is sufficient to explain the observed positive 
spectral curvature in this band. 

• At specific discrete locations in the disk, it is 
clear that the thermal contribution is not negligi- 
ble. An example is the location of a specific star 
forming region identified with a cross in Fig. [2t 
whose flat 1.5 GHz spectral index is consistent 
with thermal emission. 
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(a) X-ray over L-band 




12''42'^25' 15^ 10' 05' 00' 41 ""55= 50^ 

J2000 Riqht Ascension 



(b) L-band only 




J2000 Right Ascension 



(c) L-band over Optical 




12^^2'^-48^ 36' 12^ 00^ 4^'^48^ 36^ 24^ 

J2CfiG Rlgh: Ascenslor 



Fig. 1. — L-band (1.5 GHz) images of NGC 4631. The beam is sho wn as the filled ellip se at lower left. See 
Table [2] for map parameters (a) : Contours of soft X-ray emission (|Wang et al. I l2QQlf ) over the Robust 
weighted greyscale image, (b): Contours of the Robust 2 weighted image over a self-grey scale. Contours are 
at 9 (2cr), 13.5, 20, 30, 45, 65, 90, 150, 300, 500, 750, 1500, and 2500 x 10 /iJy beam-^ Red dashed curves 
denote the loops discussed in Sect.O (c): Contours of the uv tapered image over the Second Digitized Sky 
Survey (DSS2) blue image. Contours are at 2 (2cr), 3, 5, 8, 15, 30, 65, 120, 250, and 500 x 100 /iJy beam-^ 
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12^42'^25* 20^ IS" 10^ Q5' 00* ,41 ""55^ 50* 45* 



Right Aacanaion 



(c) C-band over Halpha 




12'^*2!^a5* 20* 15*" 10' 0!^" 4l"*5S" SQ' 



Fig. 2. — C-band (6 GHz) images of NGC 4631. The beam is shown as the filled ellips e at lower left. See 
Table [2] for map parameters, (a): Contours of CO(J=3-2) emission (Irwi n et al. ll2Qllh over the Robust 
weighted image, (b): Contours of the Robust 2 weighted image over a self-greyscale. Contours are at 11 
(2cr), 16.5, 30, 50, 80, 150, 30 0, 500, and 900 ^Jy beam'^ (c): Contours of the uv tapered image over a 
greyscale of Ha emission from lHoopes, Walterbos fc Rand I ( 19991 ). Contours are at 16 (2cr), 24, 40, 60, 100, 
200, 400, 800 and 2500 /iJy beam~^. The cross marks the location of a star formation complex discussed in 
Sect. [331 
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12''42"'24^ 18^ 12^ 06' 00' 4^54^ 48' 
J2000 Right Ascension 



Fig. 3. — Polarization maps of NGC 4631. The linear polarization is shown in greyscale with the greyscale 
range shown at right. Contours show the total intensity emission, and red vectors show the apparent angle of 
the magnetic field. (Note that the vector length could only be specified as a constant in the CASA package.) 
Maps have been corrected for ionospheric Faraday rotation, but not for Faraday rotation intrinsic to the 
galaxy itself The beam is shown as a filled ellipse at lower left, (a) L-band polarization corresponding to 
Fig. [Ud. (b) C-band polarization corresponding to Fig. [2t. 
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Fig. 4. — Map of percentage polarization at C-band superimposed with the apparent magnetic field vectors 
as shown in Fig. [Sb- The greyscale range is shown at right and the beam ellipse is shown at bottom left. 
Notice the sharp change of vector orientation at RA ^ 12^ 42"^ 12% DEC ^ 32° 32' 30''. 
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-2-1 1 




12 42 30 15 00 41 45 30 

RIGHT ASCENSION (J2000) 



Fig. 5. — L-band spectral index maps corresponding to Fig. [T^, shown to the same greyscale and spatial 
scale. The spectral index and curvature (when fitted), averaged over the displayed emission, are denoted at 
upper right. The beam is indicated as an open ellipse at lower left. The per-pixel error bar is 0.4 fSect. [216]) . 
creating the 'mottled' appearance; extreme values within a beam width of the perimeter are artifacts, (a) 
Fit to Eqn. [5] assuming no curvature {/3 = 0). (b) Fit to Eqn.O solving for both a and /3. The cross marks 
the location of a star forming complex also marked in Fig. [2t. 
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12 42 30 15 00 41 45 

RIGHT ASCENSION (J2000) 



Fig. 6. — L-band spectral index maps of Fig. [TJd and c, shown to the same greyscale (shown at top) and 
spatial scale. The spectral index and curvature, averaged over the displayed emission, are denoted at upper 
right. Contours are of the corresponding total emission. The beam is indicated as an open ellipse at lower 
left, (a) Maps corresponding to Fig. [1]d. Contours are at 9 (2cr), 20, 40, 65, 125, 250, 600, 1200, and 2400 
X 10"^ Jy beam-^ (b) Maps corresponding to Fig. [It. Contours are at 2 (2cr), 4, 8, 16, 35, 70, 150, 300 
and 600 x 10~^ Jy beam~^. 
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Fig. 7. — L-band contours of Fig. [Tt over a colour soft X-ray image of Fig. [TJi enhanced to show the X-ray - 
radio continuum association in the loop north of the galaxy's center. 
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Table 1 

NGC 4631 Observing Parameters 



Parameter 


1.5 GHz 


6.0 GHz 


Array 


C 


C 


No. antennas^ 


25 


22 


AC Central Frequency (GHz)'^ 


1.375 


5.5 


BD Central Frequency (GHz)^ 


1.625 


6.5 


Spw bandwidth (MHz) 


32 


128 


No. spws (AC plus BD)^ 


16 


16 


Total bandwidth (AC plus BD)(GHz) 


0.512 


2.048 


No. channels/spw 


64 


64 


Total no. channels 


1024 


1024 


Channel separation (MHz) 


0.500 


2.00 


Integration time (s)*^ 
Obs. Time (min)^ 


10 


10 


30 


75 


Flux calibrator^ 


3C286 


3C286 


Zero-polarization calibrator 


OQ208 


OQ208 


Phase calibrator 


J1221+2813 


J1221+2813 



^The total number of antennas for which data were acquired during 
the observing session, after problem antennas were flagged. 

^AC and BD IFs are contiguous in frequency. 

"^Number of spectral windows in AC or BD IFs. 

"^Single record measurement time. 

®On-source observing time before flagging. 

^This source was also used as the bandpass calibrator and for deter- 
mining the absolute position angle for polarization. 
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Table 2 
NGC 4631 Map Parameters 



Parameter 




1.5 GHz 






6.0 GHz 




uv weighting^ 
uv taper (kA)'^ 
Nterms'' 
No. Self-cals^ 


Briggs Rob = 
none 
3 

1 a&p 


Briggs Rob = 2 
none 
3 

1 a&p 


Briggs Rob = 2 
5 
3 

1 a&p 


Briggs Rob = 
none 
3 

1 a&p 


Briggs Rob = 2 
none 
3 

1 a&p 


Briggs Rob = 
10 
3 

1 a&p 


I images 


Figure label 
Synth, beam^ 

(", ", °) 

rms (/iJy beam 


Fig.lila) 

11.16, 10.28, -175.0 
28 


Fig.lilb) 

16.69, 14.77, 34.5 
45 


Fig.lilc) 

29.90, 26.08, 48.9 
100 


Fig.l2la) 

2.71, 2.63, -50.9 
4.8 


Fig.l2lb) 

4.07, 3.80, 27.5 
5.5 


Fig.l2lc) 

8.60, 8.52, 41.5 
8.0 








Q & U ima 


-ges^ 






Synth, beam* 

rms {fiJy beam"""^)^ 


11.01, 10.13, -176.8 
27 


12.28, 14.41, 33.7 
22 


29.35, 25.51, at 48.3 
23 


2.71, 2.63, -50.8 
4.7 


4.07 X 3.80 at 27.5 
3.5 


13.15, 11.45, 27.2 
7.0 



Note. — These parameters represent the images presented in Figs, fll and [21 
^See lBriggs I (jl995 ) for a description of Briggs weighting with various 'robust' factors. 

Scale length of Gaussian taper applied in the uv plane. 
'^Number of terms in the Taylor expansion for fitting the spectral index. 

"^Number of self-calibration iterations, where 'a&p' refers to amplitude and phase together. 
®Map labels, as shown in Figs. [Tland [2l 

^Synthesized beam major and minor axis and position angle. 
^Rms map noise before primary beam correction. 

^The cross-hands (RL, LR) had different sets of flags applied than the parallel hands (RR, LL), leading to differences in the synthesized beams and 
noise for Q and U compared to I. Cleaning also proceeded with different scales. 



